End-point of the rp process and periodic gravitational wave emission 
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The general end-point of the rp process in rapidly accreting neutron stars is believed to be a surface 
distribution of matter whose nuclear composition may depend on position. Its evolution during 
compression beyond the neutron-drip threshold density p n£ j is determined by the presence of nuclear 
formation enthalpy minima at the proton closed shells. At p n d, a sequence of weak interactions 
with capture or emission of neutron pairs rapidly transform nuclei to the most accessible proton 
closed shell. Therefore, angular asymmetries in nuclear composition present in accreted matter at 
p n d are preserved during further compression toward densities ~ 10 14 g cm -3 provided transition 
rates between closed shells are negligible. Although it has been confirmed that this condition is 
satisfied for predicted internal temperatures and for the formation enthalpy distribution used in 
this work, it would not be so if the true enthalpy differences between maxima and minima in the 
distribution were a factor of two smaller. For this reason, it does not appear possible to assert with 
any confidence that position-dependent surface composition can lead to significant angle-dependence 
of the equation of state and to potentially observable gravitational radiation. The effect of non-radial 
internal temperature gradients on angle-dependency of the equation of state is also not quantifiable. 

PACS numbers: 97.60.Jd, 97.80.Jp, 95.55.Ym, 26.60.+C 



I. INTRODUCTION 



Publication of the first data from the GEO 600 and 
LIGO interferometric gravitational wave detectors Q] is 
a notable event in the study of neutron stars. Upper lim- 
its for periodic emission were obtained, initially for PSR 
J1939+2134, and then for a further 27 isolated pulsars. 
The current rate of improvement in sensitivity makes 
the detection of signals from these and similar galactic 
sources an interesting future possibility. Gravitational 
wave emission occurs if the rotating neutron star has 
mass-quadrupole moment tensor components that are 
time-dependent in the observer frame. Magnetic stress 
is an obvious possible source of such deformation, par- 
ticularly if the internal fields are one or more orders of 
magnitude larger than those inferred at the surface. A 
brief recent survey of this topic has been given by Cut- 
ler Furthermore, it has been recognized, within the 
last ten years, that high rates of mass accretion in low- 
mass X-ray binary systems may also lead to neutron star 
deformation and gravitational wave emission [4], |5j, |fj . 

Burning of accreted hydrogen and helium in the neu- 
tron star atmosphere has been the subject of a very large 
number of published papers. Ref. gi ve re- 

sults from more recent calculations of the distributions of 
nuclear mass number A and charge Z which are formed 
by rapid proton capture (the rp process) during stable or 
unstable burning. The predicted distributions are func- 
tions of the thermodynamic conditions existing during 
the rp processes, indicating the possibility of non-trivial 
differences in atmospheric composition between, for ex- 
ample, the magnetic polar and equatorial regions of the 
star. The accretion rates considered can be extremely 
large (~ 10~ 8 M Q yr~ x ) and the consequent inward flow 
of matter is so rapid that any dependence of atmospheric 



composition on the angular coordinates 9, <fi which define 
position on the neutron-star surface must certainly be 
transferred to the liquid ocean which is present above 
the solid phase of the crust. Unless there is very strong 
convective instability it will survive there for times ~ 1 
yr before entering the solid phase at a matter density 
p ~ 10 9 g cm~ 3 . (Fig. 1 of Ref. Q gives a useful sum- 
mary of relevant orders of magnitude for surface layer 
properties.) An element of accreted matter moving in- 
ward is compressed and undergoes a sequence of nuclear 
transitions as the density p increases. Finally, there may 
be changes of phase before entry to the liquid core of the 
neutron star. Sato first recognized that the composi- 
tion of this accreted matter, as a function of p, would be 
quite different from that of primordial neutron-star mat- 
ter. Detailed calculations of the sequence of transitions 
starting from (A, Z) = (56,26) at p = 10 8 g cm~ 3 have 
been mad e by Haensel and Zdunik [12| who also repeated 
this work [Rj for the initial condition (A, Z) = (106,46) 
following publication of Ref. 0. Below the neutron- 
drip threshold p n d, matter consists of nuclei embedded 
in a relativistic electron gas. The transitions in this re- 
gion are a sequence of electron capture reactions which 
reduce Z at constant A. Thus any position-dependent 
differences in nuclear mass number present in the at- 
mosphere and surface ocean are preserved at p < p n d 
so that the equation of state may have, in principle, a 
small 9, </>-dependent component 0, 0|. Although this 
region contains only ~ 10 -5 M , the neutron-star mass- 
quadrupole tensor components generated by it are not 
necessarily negligible. Even if surface ocean composi- 
tion were uniform, a non-radial temperature gradient 
Vg^T ^ present within the solid crust changes the 
densities at which the electron capture transitions occur 
and so leads, in principle, to a small angle-dependent 
component in the equation of state 0, 0] . 
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The inner part of the solid crust at p > p n d is more 
massive (~ 10 -2 M Q ) and therefore the investigation of 
possible angle-dependency in composition and equation 
of state in this region is a very important problem. It 
has been considered recently by Ushomirsky, Cutler and 
Bildsten [§J in calculations based on the work of Sato 
[IT| and of Ref. [l^. The analysis of the evolution of 
nuclear charge at p > p n d given in the present paper is 
in substantial disagreement with Ref. [5| because these 
authors assume, following Sato 0] and Ref. 0, that 
accreted nuclei at p > p n & follow an evolutionary path of 
decreasing Z, produced by successive electron captures, 
until transition rates for pycnonuclear fusion reactions 
become appreciable at Z ks 10. Calculations of nuclear 
formation enthalpies |l4j at p > p n d, though necessar- 
ily very approximate, show that local minima exist at 
the proton closed shells with Z = 20,28,34,40,50. Our 
analysis is that the presence of these minima determines 
the evolution of accreted matter during compression from 
the neutron-drip threshold until p approaches the region 
of phase transitions to possible lower-dimensional struc- 
tures, or the liquid core of the star. The internal tem- 
perature and the nature of the formation enthalpy dis- 
tribution are the crucial factors that determine whether 
or not the system reaches weak-interaction equilibrium 
(homogeneity in Z) in the 10 6 — 10 7 yr interval of high 
accretion rates. Sec. IIA contains a brief description 
of the Z-dependence of the nuclear formation enthalpy 
which is the basis of the present paper. A system in 
which more than one closed shell is populated must even- 
tually become unstable, during compression to higher 
densities, against quantum-mechanical tunnelling of pro- 
tons between adjacent nuclei, even if weak-interaction 
transition rates are negligibly small (Sec. IIB). We also 
consider the inward movement of the accreted matter 
and show in Sec. IIIA that the stress-response of the 
solid cannot, without qualification, be described as elas- 
tic. Calculation of the quadrupole moment produced by 
the angle-dependent equation of state is therefore not a 
simple problem and for this reason the solid crust is rep- 
resented as a one-dimensional system in Sec. IIIB. A brief 
comparison is given with the one-dimensional treatment 
considered in Ref. 



Whether or not periodic signals are seen by interfer- 
ometric detectors, it is important to have some under- 
standing of the ways in which they might be generated. 
Magnetic stress is the more simple source, although ques- 
tions such as the form of the field distribution and the 
nature of proton superconductivity present themselves. 
The main conclusion of this paper is that our present 
understanding of the relevant aspects of nuclear struc- 
ture in the neutron-drip region is such that predictions 
of gravitational wave emission arising from compositional 
asymmetries cannot be made with confidence. But the 
significant factors are considered in Sec. IV. 



II. THE STRUCTURE OF THE 
NEUTRON-DRIP SOLID 

A. Shell effects above the neutron-drip threshold 

The canonical view of the solid phase has been that 
strong and weak-interaction equilibria are maintained 
during cooling so that a Coulomb bec lattice of nuclei 
is formed on solidification, neutralized by a relativistic 
electron gas. The matter density at the neutron-drip 
threshold has some dependence on the details of the equa- 
tion of state but is p n d ~ 4 — 8 x 10 11 g cm -3 . At 
P > Pnd, the neutron chemical potential exceeds its rest 
energy and the nuclei are in equilibrium with both the 
electrons and a continuous neutron medium, superfluid 
at temperatures below T ~ 10 9 K. The early and clas- 
sic paper of Negele and Vautherin |15| , involving micro- 
scopic calculations of the single-particle states for neu- 
trons and protons inside a Wigner-Seitz cell, gave the 
structure of nuclei within the neutron medium. The en- 
ergy differences per nucleon, at fixed baryon density, be- 
tween lattices with different proton closed shells, for ex- 
ample, Z = 40, 50 were found to be very small, indicat- 
ing that the standard assumption of a homogeneous bec 
lattice might be questionable. Much later, calculations 
of the formation enthalpies for monovacancies and for 
charge-impurity point-defects [l6l IT?! produced values 
which showed that the assumption could not be valid. 
In that work, a simple procedure was used to obtain 
the properties of homogeneous lattices as a basis from 
which the enthalpy changes produced by perturbations 
(the point-defects) could be calculated. Nuclei were de- 
scribed by the compressible liquid-drop model (CLDM) 
of Lattimer et al |18| with their model parameter val- 
ues. Table I contains the calculated lattice parameters 
used. The CLDM nuclear charge Z is a continuous vari- 
able, and the neutron and proton chemical potentials are 
given with reference to the neutron rest energy. The for- 
mation enthalpy for a point-defect of charge Z can be 
expressed as Hpz = C[Z — Z) 2 relative to the homoge- 
neous lattice. It is reasonable to expect that the order 
of magnitude of Hpz should be related with the mono- 
vacancy formation enthalpy Hp v through the expression 
H FZ « {ZZ- 1 - l) 2 H Fv , with C w H Fv Z- 2 . The val- 
ues of C contained in the final column are broadly con- 
sistent with this, given that H Fv w 15 MeV 0. These 
Hpz values, determined solely by Coulomb lattice and 
bulk nuclear matter properties, were small and it became 
obvious that shell corrections should not be neglected. 
Their inclusion [l4|, based on the single-particle levels 
found by Negele and Vautherin and the Strutinski 
procedure (see Ref. pjfh changes the Hfz to the for- 
mation enthalpies shown in Fig. 1. These are exclusive 
of integral multiples of the electron and neutron chemi- 
cal potentials. The quadratic dependence and 
the effect of proton closed shells are both obvious. We 
refer to Ref. for further details of the methods of 
calculation of these Hpz and of the approximations and 
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assumptions on which they are based. They are nec- 
essarily subject to some uncertainty and it is also the 
case that, for laboratory nuclei far from stability, the de- 
pendence of the spin-orbit interaction on neutron excess 
remains a problem of current interest (see Schiffer et al 
pcj)- Although the formation enthalpies given in Fig. 1 
are, as we have stressed, of uncertain quantitative value 
(which is why the simplification of approximating free en- 
ergy by enthalpy has been made), it is fairly certain that 
shell corrections producing local minima in Hpz must 
exist, even if they are smaller than those of An 
additional factor supporting this conclusion is the belief 
that the quadratic dependence of Hpz on Z — Z has 
not been seriously under- estimated. (The reason for this 
is that the monovacancy formation enthalpies obtained 
in [l6| are large when expressed in units of the melting 
temperature, Hp v ~ 35/csT m , whereas for alkali metals 
they are rs 10— 15/csT m .) Therefore, the qualitative fea- 
tures of Fig. 1, and specifically the existence of large 
potential barriers between closed-shell values of Z, are 
probably well-founded. The consequence is that weak- 
interaction equilibrium fails during cooling of primordial 
neutron-star matter |14| producing a Z— heterogeneous 
and amorphous solid. This is the primordial state of the 
solid in an isolated neutron star or in a binary system 
that has not yet evolved to a high rate of mass accretion. 



B. Elimination of heterogeneity by 
quantum-mechanical tunnelling 

The metastable state of matter described in Sec. IIA 
is not possible at matter densities approaching the phase 
transition to the core liquid. Lower-dimensional phases, 
if they exist at densities between those of the spheri- 
cal nuclear phase and liquid core |2l|, are unable to 
support point-defects. But since such phases may not 
exist, it is necessary to investigate the stability of the 
Z— heterogeneous solid for values of p larger than those 
of Table I. At finite temperatures, quantum-mechanical 
tunnelling of single protons and of pairs can occur be- 
tween any nucleus and its nearest neighbours. Consider 
two nuclei with proton closed shells Z\^- The tunnelling 
of a single proton forms particle and hole excitations, 
Zi — » Zi±l and Z 2 — > Z 2 T 1 with an associated enthalpy 
change 5Hp > which must be balanced by interac- 
tion with other degrees of freedom present in the system. 
(Pair tunnelling would allow formation of even— Z nuclei 
but with much smaller transition rates.) The nature of 
the potential barrier in this unusual system is somewhat 
uncertain, but we have assumed that it is determined by 
the difference between that proton chemical potential \i e v 
which is defined by proceeding to the limit of zero pro- 
ton number density within the neutron continuum, and 
the chemical potential pt p in the nuclear volume. The 
values of pf v given in Table I have been obtained from 
the Skyrme bulk nuclear matter pseudo-potential used by 
Lattimer et al ^| . The order of magnitude of the bar- 
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FIG. 1: Formation enthalpies for charge-Z point defects, from 
the work of [l^| , are shown for the homogeneous lattices whose 
properties are given in Table I. Integral multiples of the elec- 
tron and neutron chemical potentials are excluded. There are 
local minima at the proton closed-shells Z — 20, 28, 34, 40, 50. 



rier penetration factor, calculated for the homogeneous 
lattice of charge Z, in which h 2 k 2 = 2m p ([ip — fj, p ) and 
x = 2(r\vs ^ t n), is given in Table I. 

Scattering of neutron superfluid quasiparticles by pro- 
tons is the most important interaction with other degrees 
of freedom. A very approximate estimate of the transi- 
tion rate for the tunnelling process Z\ — » Z\ ± 1 and 
Z 2 — > Zi =p 1 can be found from the Table I barrier pen- 
etration factor. It is of the order of 
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where ft^ 1 — ksT and the parameters kp n and A„ are 
the neutron Fermi wavenumber and energy gap within 
the nuclear volume. The order of magnitude of the 
strong-interaction scattering matrix element is V ~ 10 2 
MeV fm 3 , leading to a transition rate 
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TABLE I: Homogeneous CLDM lattice properties are given for a set of matter densities p. These are the neutron continuum 
number density n„, the nuclear charge Z and radius rjv, the Wigner-Seitz radius rws, the neutron and proton chemical 
potentials p n , P and the pressure P. Also given is the CLDM proton chemical potential pp defined in the limit of zero proton 
number density within the neutron continuum. The final two columns give the proton barrier penetration factor exp( — kx) 
and the parameter C. 
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at 10 9 K. For temperatures expected either during accre- 
tion or early in the life of an isolated neutron star, and for 
values 5 Hp ~ 1 MeV that are consistent with Fig. 1 , the 
most significant exponential term in Eq. (1) is certainly 
the barrier penetration factor. Its rapid increase from 
negligibly small values is an unavoidable consequence of 
the variations of both rws and p e v — p p with increasing 
p. At values of p not much larger than the third row of 
Table I, the proton tunnelling transition rate must be- 
come so large that the metastable equilibria represented 
by the minima in Fig. 1 become short-lived and all nu- 
clear charges approach that of the most favoured closed 
shell, probably Z = 40 

Our conclusion is that the state of matter at this 
density and beyond is independent of its past history 
and that, in consequence, any direct transition from the 
spherical nuclear phase to the liquid must occur at a 
unique pressure P c and density p c . The same state- 
ment can be made of any intermediate transition to 
lower-dimensional solid phases which can support some 
anisotropic stress-tensor components [22^. This has the 
consequence that, within the Cowling approximation, 
in which the perturbation to the gravitational poten- 
tial caused by an angle-dependent solid crust composi- 
tion is neglected, the static structure of the neutron star 
is very simple. The surface separating liquid and solid 
phases is necessarily isobaric and spherical with radius 
R c , and these conditions are maintained by strong and 
electromagnetic interactions. This statement, of course, 
neglects rotation and the effect of thermal fluctuations 
which average to zero over macroscopic areas. The ex- 
tent to which it would be true for the various modes of 
oscillation of a neutron star is not considered here. 



III. EVOLUTION OF ACCRETED MATTER 
A. The inward flow of matter 

The predicted nuclear mass numbers produced in the 
rp process burning of accreted matter extend to A w 110 
00 EH • Widths and shapes of these mass num- 
ber distributions depend on thermodynamic conditions 



which are determined principally by the local rate of ac- 
cretion and its hydrogen-helium composition. Therefore, 
it is very likely that the mean mass number of nuclei 
entering the surface ocean of the star is a function of 
surface position specified by the angles 0, cj>, with the 
possibility of significant differences between the magnetic 
polar and equatorial regions. The extent to which con- 
vective motions exist in the ocean is a complex prob- 
lem but is known to depend on magnetic field strength 
and on the presence of a non-radial temperature gradi- 
ent [23, 24]. The more conservative assumption is that, in 
general, mixing is incomplete over the whole surface area 
and that some composition angle-dependence remains as 
matter enters the solid phase of the crust. This inward 
movement of individual nuclei occurs at constant A, but 
with Z decreasing owing to electron capture transitions. 
Given the width of the predicted distributions, the solid 
structure formed must be amorphous, with any order lim- 
ited to very short lengths 0] . 

However, further structural rearrangement must occur 
at a very high rate during compression because the aver- 
age nuclear cell radius decreases during the inward move- 
ment. Below the neutron-drip threshold, it is oc p -1 / 3 
and at p > p nc i, the density- variation obtained from col- 
umn 5 of Table I shows almost the same dependence. 
These variations are very much more rapid than those 
which accompany the inward radial movement of a lat- 
tice without change of form and necessarily imply that 
structural rearrangement occurs. Its effect can be seen by 
assuming that, at some instant of time, the strain com- 
ponents are e^, defined relative to suitable fixed coordi- 
nates. Given the p-dependence of rws , it follows that nu- 
clear rearrangments consequential to inward movement 
can, in principle, change the strain at a rate as large as 
&y ~ v/L p where L p ~ 10 4 cm is the scale-length for 
radial density variation and v is the bulk inward radial 
velocity of matter at density p. Although we are unaware 
of any strictly analogous terrestrial system, the high rate 
of nuclear rearrangement has some similarity with visco- 
elasticity. The solid system considered here must have 
the same elastic stress-response as an amorphous solid for 
times several orders of magnitude shorter than eijL p /v, 
but over longer times, the strain is continuously changing 
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owing to the inward movement of accreted matter. The 
important feature is that the solid, unlike a conventional 
visco-elastic structure, has no intrinsic relaxation time: 
the values of are determined by the accretion rate. 
Later in the evolution of the binary system, when the 
accretion rate has become neglible, the structure is again 
that of an amorphous solid. 

Below p ni i, the equation of state of neutron star matter 
is polytropic, P — Kp T , with an adiabatic index T very 
close to 4/3 because the pressure P is well approximated 
by the relativistic electron pressure P e . In this region, an 
incremental change in nuclear composition from reference 
mean values A and Z produces a simple change in the 
constant K. At a fixed P e , the matter density p cx A/Z 
to a high level of approximation. The adiabatic index is 
constant and therefore, 



SK S(A/Z) 
K (A/Z) ■ 



(2) 



No similar elementary result can be derived at p > p nc i 
because the neutron continuum contributes the greater 
part of the pressure. But the CLDM equation of state 
underlying Table I is approximately polytropic, with adi- 
abatic index T — 1.29, and we shall assume that compo- 
sition changes in this region also give K — > K + SK with 
no incremental change in T. 



B. Angle-dependent structure of the solid crust 

Given the considerations of Sec. IIB and III A, the ba- 
sic properties of the crust in an accreting neutron star can 
be found very simply using a polytropic equation of state 
and the Cowling approximation, in which perturbations 
of the gravitational potential are neglected. Suppose that 
the nuclear composition is angle-dependent, giving an 
equation of state with K — > K + 5K(9, <fi), but with the 
length-scale for non-radial variation of SK at least of the 
same order of magnitude as the crust depth. A three- 
dimensional elastic-deformation treatment of the crust 
is appropriate for calculations of the adiabatic modes of 
oscillation such as those of McDermott, Van Horn and 
Hansen [2{| at zero or negligible accretion rate. Nev- 
ertheless, given the many uncertainties in the formula- 
tion of the present problem, a one-dimensional (radial) 
approximation for local values of 9, <fi has the merit of 
transparency and is adequate. 

In the three-dimensional treatment, Eulerian devia- 
tions from the reference system (SK = 0) satisfy the 
full static equilibrium condition (or quasi-static in view 
of the accretion flow) 
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+ 8pg l = 0, 



in which the deviation in the stress tensor, 



So, 



SPSij, 



(3) 



(4) 



is the sum of the deviation in the shear components of 
the elastic stress tensor, the Maxwell tensor, and the 
term — SP which includes the deviations in all other 
isotropic stress tensor components of the system. Let 
us assume that SK ^ within a pressure interval < 
P < P a , where P a < P c . (It is possible that P a may 
be slightly composition and therefore angle-dependent, 
but this would have an effect one or more orders of mag- 
nitude smaller than those of first order in SK.) In the 
interval P a < P < P c where matter has reached equilib- 
rium, either through high rates of proton tunnelling or 
a transition to a lower-dimensional structure, the angle- 
dependence vanishes so that SK = 0. Eq. (3) satisfy the 
boundary conditions Sa rr = at pressures P = 0,P a , 
though possibly with a small error if P a < P c . The den- 
sity deviation, integrated from R a to the surface of the 
star is then 



Spdr 



R d5o r 



'1r JR a dXj 



-dr. 



(5) 



In a true one-dimensional approximation, the integrated 
density deviation therefore vanishes. It is worth compar- 
ing Eq. (5) with Eq. (51) and (52) of Ref. @ which 
satisfy identical boundary conditions. But these authors 
find, by reference to their three-dimensional solutions, 
that | Sa rr | is at least an order of magnitude larger than 
| Sa r j | throughout most of the interval of integration. 
They therefore conclude that a one-dimensional estimate 
of a quadrupole moment which they express, in Eq. (52) 
of their paper, as 



Q22 = R 4 / Spdr 



(6) 



is extremely poor owing to its independence from the 
major stress deviation Sa rr . It is obvious that this is 
correct. (In the one-dimensional approximation which 
will be adopted here, the moment so defined would be 
Q22 — 0.) These authors observe that changes in the 
radial distribution of density produced by SK ^ are a 
far larger contributor to quadrupole moments than the 
non-radial displacements which also arise. Although the 
latter can be large (~ 10 4 cm), they are smaller than 
the typical scale length for non-radial variation of SK 
and so are much less significant. For a thin crust, it 
must be the case that any calculation which gives the 
correct radial distribution of density deviations also gives 
a satisfactory approximation to the quadrupole moment. 
A one-dimensional solution of the form given here in Eq. 
(7) and (8) is capable of this. But it is necessary to use 
the correct quadrupole moment definition given by Eq. 
(1) of Ref. || instead of Eq. (6). 

Solutions of Eq. (3) are easily obtained in the one- 
dimensional approximation for the case in which SK is 
independent of depth within the pressure interval < 
P < P a . We replace Sa rr by — SP and find that the 
incremental changes from the reference values of pressure 
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and density are, 




r-i 



r-i 



KT(T - 1 



(7) 



The radial displacement of a point fixed in Lagrangian 
coordinates is, 



5r W 



SK 



g a (r 



V) (Pa ~ P T1 ) 



(8) 



where g a is the magnitude of the gravitational acceler- 
ation g at pressure P a . There is a density discontinu- 
ity — p a SK J KT at P a . The radius of the stellar surface 
(P = 0) is then angle-dependent and the case SK > 0, 
a slightly more stiff equation of state, produces a local 
increase in radius. 

Eq. (4) makes obvious the fact that the term SP 
produced by an angle-dependent SK ^ has an effect 
analogous with the isotropic components of the Maxwell 
tensor. The change of phase to the liquid occurs at a 
fixed pressure P c and density p c . In the Cowling ap- 
proximation, and neglecting the effect of rotation, the 
radius R c is angle-independent, but its magnitude re- 
duces very slowly with time for most equations of state 
during accretion as the total mass of the star increases. 
Solution of the three-dimensional Eq. (3) has not been 
attempted. It would not be simply a problem in elastic- 
ity because is time-dependent for non-zero accretion 
rates as we have described in Sec. IIIA. From Eq. (3) 
and (4), and with neglect of the Maxwell tensor, we can 
see that the non-radial components of the shear stress de- 
viation Scrfj are of the same order of magnitude as SP. 
Hence the maximum SP for which the solid will respond 
quasi-elastically is of the order of 10~ 2 p, for an amor- 
phous structure, where p is here the shear modulus. For 
greater values of SP, the behaviour of the system will 
approach that of a fluid. The very complex problem rep- 
resented by this condition is not considered here. In the 
neutron-drip region, the pressure and shear modulus have 
almost identical p-dependences and so are approximately 
related by p « 6 x 1CT 3 P. Therefore, the maximum SK 
for which the solid behaves quasi-elastically is given by 
SK/K » \Q- 2 ~p/P w 6 x 1CT 5 . We can see that the 
inward movement of matter during accretion is basically 
radial but, given the non-radial gradient in SP, it cannot 
be exactly so. 



C. Nuclear transitions above the neutron-drip 
threshold 



The neutron-drip region at p > p ne [ contains almost all 
the mass of the crust. Some information about the dis- 
tribution of A and Z for matter moving into this region 
has been provided by Haensel and Zdunik (see Fig. 1 



of 0). Successive electron captures in the outer crust, 
for the two mass number examples A = 56, 106 consid- 
ered by these authors, lead to the formation of nuclei 
with charges Z — 18, 32, respectively, as the inward mov- 
ing accreted matter is compressed to p n d- The assump- 
tion which appears to have been made by Sato [Tll | and 
by later authors 0] is that, with further compression 
inside the neutron-drip region, electron capture transi- 
tions continue at constant A until the transition rates 
for pycnonuclear fusion reactions become appreciable at 
Z w 10. Comparison with calculations of the neutron- 
drip state, for example, those of Negele and Vautherin 
|15| , shows that the charges found by Haensel and Zdunik 
(Z = 18, 32) lie well below the equilibrium Z. Therefore, 
unless some special case for the contrary can be made, 
the nuclear transitions which occur must follow a path of 
increasing Z. They are, principally, the successive cap- 
ture of neutron pairs followed by electron emission. For 
p n > 0, nuclei can be viewed as bound states of pro- 
tons embedded in a neutron continuum. Thus there is 
no barrier to increasing nuclear A and Z values by these 
transitions. Pycnonuclear fusion rates are negligible ex- 
cept possibly at the low-.Z end of the nuclear number 
density distribution. 

The work of Ref. ^3] shows that the inward moving 
accreted matter may have a fairly wide distribution of Z 
at p n d , but with an rms value Z a which is probably much 
smaller than either the equilibrium values calculated by 
Negele and Vautherin or those given in Table I. In or- 
der to see how the accreted matter changes as it is further 
compressed to p > p n d it is convenient to use, with some 
modifications, the CLDM procedure which is the basis 
for Table I. Strong and electromagnetic-equilibrium con- 
straints are retained, but the condition Spi = 0, where 



Spi = p n - Ppi - Pe ~ {m p - m n )c 2 



dfc 
dZ 



(9) 



and f c is the Coulomb energy of the Wigner-Seitz cell 
[l8j . is removed so that the nuclear charge Z can be 
treated as a constant parameter. For the modified CLDM 
procedure, the distribution of Z-values is approximated 
by a binary system (i = 1,2). The proper- frame neutron 
chemical potential p n is related with the time-like com- 
ponent of the metric and so, to a good approximation, 
can be assumed a fixed quantity at any point. However, 
the electron chemical potential p e , defined here as the 
Fermi energy, necessarily has a common value for the bi- 
nary system which has to be determined locally. This 
has been achieved by modifying the constraint 0] re- 
lating the nuclear surface thermodynamic potential den- 
sity, defined here as cfj, with the Coulomb energy of the 
Wigner-Seitz cell. (In the binary system and as a matter 
of convenience, we refer to the electrically neutral sphere, 
with nucleus Zi at the origin, as the Wigner-Seitz cell.) 
Clearly this cannot be satisfied simultaneously for the 
binary mixture i — 1, 2 and so we have adopted the form 



a,i {2Trr 2 Ni d l - f cl ) = 0, 



(10) 
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for this constraint, where r^i is the nuclear radius and 
a,i the number density fraction . These are the two mod- 
ifications to the CLDM procedure that have been made. 
The two remaining constraints in the CLDM procedure, 
those of neutron chemical potential and nuclear pressure 
equilibria are satisfied by both nuclei of the binary sys- 
tem. 

Although the modification given by Eq. (10) is a very 
elementary approximation for the binary system, the 
most important strong and electromagnetic-interaction 
constraints have been retained. Calculations of the equi- 
librium, using this CLDM procedure, confirm that the 
basic properties of accreted matter change in an intu- 
itively obvious way as it is compressed to p > p n d- For 
the case in which the binary system mean square nu- 
clear charge J2 a i^f = where Z is the equilibrium 
charge given in Table I (for the lowest matter density, 
and with Z\ = 28 and Z2 = 40), we find that the com- 
puted values of Spi are very small and that there are 
negligible changes in P, p and chemical potential per 
baryon. This is consistent with the discussion of for- 
mation enthalpy differences given in Sec. IIA. But the 
mean square nuclear charge Z 2 of accreted matter at p n d 
is smaller than Z 2 . Therefore, we have computed 

binary model cases in which Z\ = Z and Z 2 = 20 and 
find that the chemical potential imbalance defined by Eq. 
(9) is an approximately linear function of a 2 , and is given 
by Sp2 ~ 0.6 + 27a2 Mev. The implication is that weak- 
interaction transitions rapidly reduce the value of a2 at 
P > Pnd and change the system toward the equilibrium 
state of Table I. Ultimately, for small values of Z—Z a and 
therefore of 02, the nuclei of charge Z2 can be regarded 
as impurities in the Z system so that the formation en- 
thalpy distribution of Fig. 1 becomes valid and the speed 
and efficiency of this process are determined by the tem- 
perature and by the formation enthalpy barriers between 
the closed shells at Z = 20, 28, 34 described in Sec. IIA. 

Nuclear reactions caused by compression of the ac- 
creted matter are a volume source of heat. Ref. p| 
gives temperature distributions for an accretion rate 
of 10~ 8 M Q yr _1 calculated using various assumptions 
about neutron superfluidity They are fairly slowly vary- 
ing functions of p with maxima of 7 — 8 x 10 s K at p n d- 
(These high values show that a classical neutron gas must 
be present in a finite region with p n < m n c 2 .) For this 
accretion rate, the inward radial velocity is 6 x 10 -8 cm 
s -1 , at the threshold p nd = 6 x 10 11 g cm -3 0. Initially, 
for large imbalances 5pi > 0, weak-interaction transitions 
have no formation enthalpy barrier and are rapid. But 
as the value of Z a increases toward Z and the forma- 
tion enthalpy distribution of Fig. 1 becomes valid, the 
question arises of the extent to which thermal excitation 
of the weak transitions allows movement between proton 
closed-shells. Low enough formation enthalpy barriers 
would allow all nuclei to reach a unique proton closed- 
shell well within the transit time of accreted matter in 
the crust. 

This is an important question because the modified 



CLDM procedure described above shows that a nuclear 
charge imbalance, of the form Z a < Z, has a significant 
effect on the equation of state. This should not be un- 
expected because both the nuclear pressure equilibrium 
constraint and Eq. (10) depend quadratically on Z. For 
the binary system charges Z\ = Z and Z2 = 20, at the 
lowest density of Table I, the incremental changes SP and 
Sp with respect to the system with mean charge Z give 
an incremental change in the equation of state, 



SK SP Sp 

_ = _-r- = 0.15a 2 , 



(11) 



for an assumed constant V = 4/3. In terms of 
quadrupole-moment generation, this would be a po- 
tentially enormous effect if its presence were angle- 
dependent. The problem becomes one of estimating the 
interval of p over which a-i is non-negligible. 

Assuming the formation enthalpies given in Fig. 1 
(lowest density), the rates for Z = 23 — > 24 and Z = 
25 — > 26 have been calculated as functions of temper- 
ature. The approximations made in calculating weak- 
interaction transition rates are precisely those described 
in previous papers [l4l Il7| . In the temperature region 
of interest, the time-constant for the depletion of a small 
value of (Z2, with Z2 = 20, is r ~ 3.3exp(51/Tg) s, where 
Tg is the temperature in units of 10 9 K. The times, for 
the temperatures predicted at p n d in Ref. [5j , are so long 
that there is no question that significant non-zero 02 can 
survive the ~ 10 s interval of rapid accretion. But we 
have previously stressed the many sources of uncertainty 
in the details of Fig. 1. The existence of these problems 
is unfortunate because halving the formation enthalpy 
barrier (roughly equivalent to halving the exponent in 
the time-constant expression) reduces the depletion time 
constant to values within one or two orders of magnitude 
of the accretion time-interval. 

In summary, the evolution of accreted matter in the 
neutron-drip region is as follows. A nucleus of charge Z 
at p n d undergoes rapid weak-interaction transitions, with 
capture or emission of neutron pairs, toward the most 
accessible proton closed shell configuration. The calcu- 
lations of burning 0, H, 0, 03 and of evolution below 
the neutron-drip threshold 13] indicate Z-distributions 
which would evolve by transitions toward those at Z = 
20, 28, 34 and possibly 40. Their relative populations re- 
flect the width and form of the distribution of A formed 
in the atmosphere. But the mean square charge Z 2 at 
p nd is appreciably smaller than Z 2 so that, as our CLDM 
binary system calculations confirm, there will be further 
rapid transitions toward one or more closed shells near 
Z until the populations of the other closed shells become 
small, of the order of aj ~ 10~ 2 . At this stage, the system 
may be metastable or, depending on temperature and on 
formation enthalpy barrier height, there may be further 
transitions to a unique closed-shell Z. As we have em- 
phasized in the previous paragraph, it is unfortunately 
the case that we cannot decide between these two pos- 
sibilities with any confidence. The sequence of forma- 
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tion enthalpy distributions shown in Fig. 1 indicate that 
compression to p > p n d may itself change the closed- 
shell populations. Thus the Z = 28 minimum becomes 
less pronounced in the vicinity of 3.7 x 10 13 g cm -3 and 
it is possible that its population is transferred firstly to 
Z = 34 and then, near 8.8 x 10 13 g cm" 3 , to Z = 40. But 
it must be re-emphasized that the specific values of Hpz 
in Fig. 1 should not be taken too seriously 0] and that 
although they show the general way in which the system 
evolves, any details that have been presented for illustra- 
tive purposes are not necessarily reliable. These changes 
under compression continue until the density reaches p a , 
where the transition rates for quantum-mechanical pro- 
ton tunnelling become so large that any metastability 
disappears. 

Non-radial temperature gradients, if present in the 
solid crust, are an independent source of composi- 
tion angle-dependency on any isobaric surface 0, |^] at 
P < Pnd- This arises, in principle, because the weak- 
interaction transition rates which change Z as p increases 
are temperature-dependent. Above p n d, the initial pro- 
gression of a nucleus toward a particular proton closed 
shell is rapid because there is no formation enthalpy bar- 
rier. Thus temperature is unlikely to be a significant fac- 
tor. But sequences of transitions between adjacent closed 
shells are, as we have described, extremely temperature- 
dependent and it is not possible to predict these rates 
with any confidence. 

IV. CONCLUSIONS 

Unfortunately, the main conclusion reached in this 
paper is that it is not possible to calculate the angle- 
dependence of the equation of state (the incremen- 
tal function SK(9, </>)) arising from either an angle- 
dependent atmospheric composition or non-radial tem- 
perature gradients in an accreting neutron star. We have 
shown that such a calculation depends on nuclear struc- 
ture properties that are very unlikely to be established 
with any degree of certainty in the immediate future. 
But this conclusion, though unfortunately negative, is 
not without value in relation to future developments in 
gravitational wave detection. 

The hypothesis that the limiting rotation frequencies of 
neutron stars in low-mass X-ray binary (LMXB) systems 
are determined by gravitational radiation 0, [B| specifies 
a fairly compact range of values for the mass-quadrupolc 
tensor component. Integration of either Eq. (5) or (7) 
shows that an increment SK ^ produces no change in 
the total mass within the interval < P < P a though 
there is a density discontinuity Sp a at P a . Therefore, the 
resultant mass-quadrupole tensor component, as defined 
by Eq. (1) of Ref. |f| with Spi m replaced by Sp^\ arises 
solely from the mass rearrangement given by Eq. (7) and 
is 

Q22 « l^PaRl 8 ^, (12) 



where h is the depth of the SK ^ layers. If these 
are confined to p < p n d, for layers such that h 2 p a « 
10 19 g cm -1 , the predicted mass-quadrupole tensor is 
several orders of magnitude smaller than the specified 
value, Q22 = 3.5 x 10 37 g cm 2 [5j. In the neutron-drip 
region, values of h 2 p a ~ 10 24 g cm -1 are possible, for 
which the specified Q22 would be given by SK k 2 x 
10~ 4 K. It is also interesting that these SK are not much 
different from the approximate value corresponding with 
the quasi-elastic limit noted at the end of Sec. IIIB. This 
is not inconsistent with the conclusions of Ref. p| . With 
reference to Eq. (11), we can see that only a very small 
charge imbalance, if suitably angle-dependent, is needed 
to produce these SK at a given density. But, as we have 
emphasized, it is not possible to predict whether or not 
the metastable population of proton closed shells can be 
sufficiently long-lived to maintain this. 

The above conclusions are for angle-dependent com- 
position asymmetries formed in the atmosphere of the 
star and make no reference to composition asymmetries 
produced by non-radial internal temperature gradients 
PI 0] which may be present. Those formed at p < p n d 
have too small an effect to be of significance but at 
P > Pnd, there is a possibility that the conclusions of 
Sec. IIIC may require modification. The value of Z a at 
p n d may have a small temperature-dependence, but this 
will not affect the population of proton closed shells ini- 
tially formed. The major effect of temperature asymme- 
try is in changing the transition rates from these closed 
shells to the closed-shell Z of complete weak-interaction 
equilibrium. As emphasized in Sec. IIIC, these rates 
are exponentially-dependent on formation enthalpy dif- 
ferences that are not well-known. Although they are 
probably many orders of magnitude too small or large 
for conceivable temperature gradients to have any effect, 
there is a non-quantifiable though small probability that 
transition rate orders of magnitude may allow the forma- 
tion of angle-dependent proton closed-shell populations 
within an appreciable interval of p. 

It is worth considering briefly the extent to which in- 
trinsic temperature-dependence of the equation of state 
can lead to asymmetry in the presence of non-radial tem- 
perature gradients. There appears to be no published 
work on this problem for the case of the solid phase of 
neutron star matter. Calculations at high temperatures 
|l8j | assume a normal neutron continuum, also the pres- 
ence of classical (translational) degrees of freedom, and 
so cannot be extrapolated to T < T m . In this region, 
the most important source is likely to be the tempera- 
ture dependence of the electron partial pressure P e . For 
the order of magnitude of temperature variation assumed 
in Ref. 5] (ST ~ 0.05T at T w 4 x 10 s K), the resul- 
tant increment of SK ~4x 10~ S K at a typical chemical 
potential p, e = 60 MeV is too small to be of significance. 

The origin of the differences between this paper and 
Ref. have been described in Sec. III. The evolutionary 
path of accreted matter at p > p n d is the major area of 
disagreement. Our analysis is that it is determined by the 



9 



proton closed-shell structure of nuclei. The disagreement 
concerning the treatment of the solid crust is much less 
important and in Sec. IIIB we have shown how it arises. 

A three-dimensional calculation of the structure for 
an angle-dependent equation of state forms a large part 
of Ref. 0]. These authors assume an elastic response 
and the analysis given draws on that of McDermott, Van 
Horn and Hansen [25|. But we have shown in Sec. IIIA 
that the elastic-response assumption is not valid during 
high accretion-rate intervals which necessarily produce 
nuclear rearrangement within the solid structure. In Sec. 
IIB and IIIB, we have shown that the basic boundary 
condition in the accretion problem is that the pressure P c 
at the crust-core boundary is a constant. In the Cowling 
approximation, and with neglect of the effects of rota- 
tion and thermal fluctuations, the boundary surface is a 
sphere whose radius R c , for most equations of state, de- 
creases slowly as the mass of the star increases. There is 
no impediment to the maintenance of this condition by 
the rapid transfer of matter between the solid and liquid 



phases through electromagnetic and strong interactions, 
and there is no immediate local dependence on weak- 
interactions transition rates. This condition leads to the 
very simple one-dimensional solutions Eq. (7) and (8) of 
Sec. IIIB which appear hydrostatic in nature. Nonethe- 
less, the stability of the system obviously depends on the 
shear modulus of the solid. 

As a source of gravitational radiation, it appears from 
our analysis that the deformation produced by angle- 
dependence of the rp process end-point is not quantifiable 
and may be no more important than that derived from 
the Maxwell tensor in the solid crust. The two sources are 
not distinguishable unless there is some a priori informa- 
tion about the internal field. The uncertainties exposed 
by our analysis seem unlikely to be removed in the near 
future and indicate that it will not be possible to decide 
with any confidence whether angle-dependent composi- 
tion or the magnetic structure and superconductivity of 
the core |2( is the more probable origin of any periodic 
signals seen in future gravitational wave experiments Q . 
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